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Abstract 


We implement an MCC solver called SSMCC, including MRV, conflict- 
directed or weighted (WTD) and failure rate-based (FRB) heuristics. We 
compare the performance of the algorithms on a wide variety of prob- 
lems. We also compare these new algorithms with Knuth’s existing DLX3 
solver. We conclude that the new SSMCC solver with the FRB heuristic 
is performing best in most cases. 


1 Dancing links and dancing cells 


Donald Knuth’s dancing links algorithm is a well-known method for solving 
exact cover problems (XCC) [1]. An implementation is available on his website 
as DLX2 [3]. He has also extended his method to allow multiplicities to be 
specified for each item. In this case, we have an MCC problem. The program 
for solving these problems is DLX3 [4]. 

Knuth’s first programs implement dancing links using linked lists. Recently 
he has also made available variations of his XCC-solver based on sparse set 
representations instead of linked lists. One of these programs is based on binary 
branching and is called SSXCC-BINARY [6]. 

In addition, he has provided two change files for SSXCC-BINARY, one that 
implements a ’conflict-directed’ or weighted heuristic (WTD) [7]. The other 
implments a failure rate based heuristic (FRB) [8]. 

The sparse set based algorithms and the heuristics are explained in [2]. 

What was missing until now was a DLX3-equivalent based on sparse set 
representations. The purpose of this note is to test a new solver, SSMCC, that 
implements the algorithm of DLX3 using sparse set representations instead of 
dancing links. This solver was implemented by Filip Stappers, borrowing heavily 
from the existing implementation by Donald Knuth. 


2 Code 


All code is available on http://archive.org/details/ssmcc. 


ssmcc.w is the main program. This program is derived from Knuth’s ssxcc- 
binary.w, but mixes in the ideas from his dlx3.w program to handle multiplici- 
ties. 

ssmcc-wtd.ch and ssmcc-frb.ch are change files for creating the versions that 
use the WTD and FRB-heuristic. 

Finally, for each of the above, there is a change file for using the sharp/nonsharp 
preference heuristic. 

There is also included a folder mcc-benchmarks, with the dlx-files for the 
benchmarks used below to compare the programs. Many of these files were 
generously made available by Donald Knuth. 


3 Comparison of the algorithms 


In the tabel below, the run times for the various algorithms in mems is shown 
and the size of the search tree (in meganodes). Some entries show question marks 
for the timings, in that case the program was running for several teramem but 
not seeming to be going anywhere. 


options, items, DLX3 MRV WTD FRB 
solutions) 
A 811, 202+0, 60568) 59.0Gmem 133.9Gmem 24,4Gmem 68.4Gmem 
11.7Mnodes 30.0Mnodes 8.8Mnodes 14.4Mnodes 
Cc 64, 65+0, 4860) 4.2Gmem 11.2Gmem 11.3Gmem 11.0Gmem 
7.3Mnodes 2.5Mnodes 2.4Mnodes 2.3Mnodes 
I 1740, 280+400, 8) 1.5Tmem stopped after 7Tmem | 657.5Gmems 376.83Gmem 
504Mnodes 181.3Mnodes | 109.2Mnodes 
M 1504, 88+102, 696) | 200.0Gmem 158.7Gmem 239.3Gmem 105.4Gmem 
33.1Mnodes 24.5Mnodes 45.0Mnodes 17.9Mnodes 
P 2436, 1730+0, 112) 451Gmem 527.0Gmem 1.03Tmems 418.5Gmem 
57.3Mnodes 75Mnodes 125Mnodes 55.0Mnodes 
Q* 3940, 65+126, 512) 284Gmem 2?.2TM 138.3Gmem 138.3Gmem 
19.4Mnodes ??Mnodes 18Mnodes 18Mnodes 
Rt* 13052, 364-46, 6) 28.1Gmem 136,9Gmem 20.6Gmem 17.4Gmem 
1.3Mnodes 17.3MNodes 1.4Mnodes 0.1MNodes 
Vv 1335, 15+61, 720) 519.4Gmem 29.3Gmem 10.2Gmem 7.5Gmem 
216Mnodes 9.0Mnodes 3.0Mnodes 2.2Mnodes 
W 2071, 447+0, 0) 2.4Tmem 62.1Gmem 3.2Tmem 7.9Gmem 
199Mnodes 6.2Mnodes 291 Mnodes 0.6Mnodes 
x* 576, 115+128, 4) 550Gmem 418.7Gmem 180Gmem 476.1Gmems 
530Mnodes 426Mnodes 148Mnodes 468Mnodes 


Most of the test problems are based on problems described in [1]. 

Problem A is the authentic partridge puzzle with n = 6 (7.2.2.1-155). Prob- 
lem C asks for 5 queens that control all squares on a chessboard (7.2.2.1, page 
93-94). Problem I is Ian Tullis’ problem from 2022, described in [9]. Problem 
M is finding all motley dissections of a 6x12 rectangle (7.2.2.1-369). Problem 
P is packing 28 3x4x5 cubes in a 28x28x28 box (7.2.2.1-350). Problem Q is 
finding all ways to fit 5 Qg in a Qs (7.2.2.1-162). Problem R is the word rect- 
angle problem described in 7.2.2.1, pages 94-95. Problem V was found in the 
knuth-dlx folder as pentominoes-2x2x3x5-noV.dlx. Problem W is Wainwright’s 
classical partridge puzzle with n = 6 (7.2.2.1-157). Problem X asks to find all 


Winner 
WTD 
DLX3 

FRB 
FRB 
FRB 
WTD/FRB 
FRB 
FRB 
FRB 


WTD 


8x8 crossword diagrams that contain exactly 12 3-letter words, 12 4-letter words 
and 4 5-letter words (7.2.2.1-111). 

Problems with a * use the sharp-preference heuristic. Problems with a ** 
use the nonsharp-preference heuristic. 

The results show FRB as the fastest algorithm in most cases, with WTD 
and DLX8 also being best sometimes. 


4 Further tests with SSMCC-FRB 


The performance of SSMCC-FRB on problem W is surprisingly good. Alas, 
trying it on the partridge puzzle with n=7 doesn’t have the same spectacular 
short solve time. Here the program takes 1.55Tmem/57.1Mnodes. 

The pentacubes-Ttower26.dlx problem (answer 7.2.2.1-343, page 512) is shown 
to have no solutions after 2.5Tmem by SSMCC-FRB (vs 5.6Tmem for SSMCC- 
WTD and 7.6Tmem for DLX3 reported on page 512). 

Exercise 7.2.2.1-300 contains several teramem-sized problems. Curiously, 
DLX3 benefits from using the nonsharp preference heuristic here, but the per- 
formance of SSMCC-FRB is worse with it. Both diagonals empty (pentominoes- 
10x10-nodiag.dlx, 4069 options, 133+0 items,31 solutions): 1.45Gmem for SSMCC- 
FRB (vs 1.87Gmem for DLX3 with the nonsharp preference heuristic). 

Both diagonals filled (pentominoes-10x10.dlx, 4068 options, 132+0 items, 
5347 solutions): 6.67Tmem for SSMCC-FRB (vs 4.6Tmem for DLX3 with the 
nonsharp preference heuristic). DLX3 comes out ahead! 

Empty spaces symmetric around the diagonals (pentominoes-10x10-diag.dlx, 
3959 options, 132+0 items, 1094 solutions): 1.58Tmem for SSMCC-FRB (vs 
19.2Tmem for DLX3). 

Central symmetry (pentominoes-10x10-180.d1x,4038 options, 132-0 items,98 
solutions): 190.4Gmem for SSMCC-FRB (vs 3.2Tmem for DLX3). 

Quadrants (pentominoes-10x10-quadrants.dlx): 13.9Tmem (13.7 with pre- 
processing) for SSMCC-FRB (vs 23Tmem for DLX3). 

Answer 7.2.2.1-346(e) asks to prove r(6, 6,6) < 15. This can be done using 2 
dix files (tripod666a.dlx and tripod666b.dlx). This results in (81.4+-47.7)Gmem 
for SSMCC-FRB vs (95.3+-47.6)Gmem for SSMCC-WTD vs (163.0+89.5)Gmem 
for DLX3. 

Exercise 7.2.2.1-369 mentions the motley dissections for 7xn are not fully 
explored. I tried finding them for n=(8, 9, 10, 11, 12, 13) and found (2023902, 
6157968, 11669721, 15716204, 16306460, 12659946) solutions in (100G, 485G, 
1.9T, 6.3T, 16.6T, 41.3T)mem using SSMCC-FRB. If the size keeps going 
roughly times 3 for each increasing n, the 7x17 case will take several thousand 
Tmem! I guess that case will not be solved on my little laptop. 


5 Conclusion 


We have succeeded in writing an MCC solver based on sparse set representations. 


SSMCC-FRB is the quickest solver on most problems, often being quite a bit 
faster than DLX3. If you have an MCC problem to solve, it seems the default 
choice to tackle the problem. 


6 Addendum 


After this note was written, Donald Knuth decided to discuss the SSMCC al- 
gorithm in The Art Of Computer Programming. He has made several improve- 
ments and bugfixes to the original program. I suggest to use his version of the 
program from now on [5] . I will leave the initial version of the program here 
just for reference. 

I also discovered that my DLX3 tests were done with the DLX3-SHARP 
program, using the sharp preference heuristic, instead of the normal DLX3. 
This caused incorrect runtimes to be reported. 

After Knuth’s modifications, SSMCC-FRB is not the clear winner anymore 
(See 7.2.2.3-(134) in [2]). 

In the meantime, I have also created an MCC-version of DLX-PRE, to pre- 
process dlx files. The code is available at https://archive.org/details/ 
ssmcc/d1x3-pre. 
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